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ABSTRACT 

From the study of X-ray light curve and color-color diagram of the low mass X-ray 
binary GRS 1915-1-105, observed by on board proportional counter array (PCA) of 
Rossi X-ray Timing Explorer (RXTE), we discover a new class of variability, which 
we name e class. We have studied observations between MJD 51200 and 51450. The 
class shows unusual periodic-like variation in count rate during rise time of two x-ray 
bursts. The class take place when the source is in radio quiet state. The huge expansion 
in color-intensity diagram indicates the class to be an adjusting stage of increasing 
accretion rate. Spectral analysis shows that during lower count rate, the spectrum 
is hard power-law dominating, indicating similarity towards hard intermediate state, 
and during higher count rate, the spectrum is thermal disk blackbody component 
dominating, indicating similarity towards high soft state. Hence, this class is important 
in understanding the way of state transition leads to change in accretion rate. No 
signature of any low frequency quasi periodic oscillation was seen in this class. We 
also find that when the class was showing higher counts, the average RMS amplitude 
is significantly high for high energy band (14-60 keV) compared to low energy band 
(2-8 keV). 
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1 INTRODUCTION 

GRS 19154-105 is one of the well observed low mass X-ray 
binary and was discovered by WATCH all-sky monitor o n 
board Granat satellite in 1992 (jCastro-Tirado et al.iri992h . 
The source is one of the brightest X-ray source in the sky 
and shows superluminal ejections (|Mirable and Rodrigue j 
1 19941 ) when the velocity of the radio jet is more than 
0.9c. GRS 1915-H05 is a b lack hole with 14.0 ± 4.4M0 
l|Harlaftis and Greineij l2004i ) which is consisting of a K- 
M III binary star of O.SM© in a wide 33.5 day orbit 
Greiner et"aL"200l'). The source is situated ~ 12 kpc far 
Mirable and Rodriguez 1994; Fender ot al. 1999). The GRS 



191 5-1-105 shows c onside rable variability in its light curve 
and iBelloni et al.l (|2000l ) showed that different patterns in 
light curve and hardness ratio of the source is repeating and 
classified the light curve of the source in 12 separate classes, 
namely a, /3, 7, S, 6, k, A, fi, v p, (j) and X- Later discovery of 
another two classes of variabili ty were reported l|Naik et al.l 
l2002l : iHannikainen et al.l [ioosi ) which are not so common 
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like other 12 classes. iBelloni et al.l l|2000t ) also phenomeno- 
logically classified X-ray states of GRS 1915-1-105 into three 
basic states, two softer states namely A, B, when the full 
accretion disk is observable and one harder state C when 
the innermost part of the accretion disk is not observable. 
The source spends most of the time in relatively hard state 
C. Here we are reporting discovery of one new class from the 
study of light curve and color-color diagram of the source, 
observed by Rossi X-ray Timing Explorer (RXTE). We are 
calling the new class as e class. 



2 OBSERVATIONS AND DATA ANALYSIS 

RXTE is regularly observing GRS 19 15-H05 from 1995 
(|Greiner et al.lll996l : IBelloni et al.ll2000l ). We have analyzed 
RXTE proportional counter array (PCA) data starting from 
MJD 51200 to MJD 51450 which includes total 97 obser- 
vations. The data were reduced using LHEASOFT version 
6.6.3. We used the standard convention of good time inter- 
val (exclusion time for the South Atlantic Anomaly (SAA) 
of 30 minutes, offset pointing < 0.02° and elevation above 
earths limb> 10°) for our observation. We first extracted 
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Figure 1. {Upper panel:) RXTE All Sky Monitor Liglit curve, 
at tlie time of the occurrence of this class. {Lower panel:) Green 
Bank Interferometer (GBI) radio flux density of the source in the 
same time. In both of these figures the time of new class is marked 
in blue color. 



rection has been applied to those spectrum files, then re- 
sponse matrix file is created from source spectrum using hx- 
trsp. Three PCA and three HEXTE files are fed to XSPEC 
for fitting. After background subtraction, we found that the 
HEXTE background subtracted data is only about 13 per- 
cent of un-subtracted spectrum and highest HEXTE count 
rate is less than 1% of highest PCA count rate found in 
spectrum. Since the count rate in HEXTE data is very low, 
errors are very high and so for spectral fitting, we did not use 
HEXTE data. Since the PCA count rate is low towards very 
low (< 2keV) and very high energy end (>30 keV), we have 
also ignored those energy for spectrum fitting and effective 
energy range of 3 — 30 keV was used for fitting. 1% system- 
atic error has been included. We have separately fitted data 
for type I and type II. We have chosen standard model of 
black hole consist ing of p owerlaw, multi-colored disk black- 
body component (iKubot a ct al. 1998) and a small gaussian 
to take care of 6.4 keV Fe line in accretion disk. For type II, 
a single power law and diskbb was needed to fit the spectrum 
but for type I, a broken power law was needed along with 
disskbb. These additive components which are sole proper- 
ties of the disk has been modified by interstellar absorption 
with column density. For two-component model fitting, the 
value of Nh is found to be 6 to 8 xlO~^^ with 95% con- 
fidence level. Hence, we have fixed the column density to 
7 X IQ^^ cm~^ for uniformity. 



light curves in three different energy bands to produce hard- 
ness ratio, color-color diagrams and hardness intensity dia- 
grams using cleaned 1-s resolution data. We define 3 bands 
as 2.0-5.8 keV (X), 6.0-13.6 keV (Y) and 14.0-60.0 keV (Z) 
and colors are defined as HRl — Y/X and HR2 = Z/X. 
For extracting light curve of different bands, we used same 
number of PCUs for all the observations. To probe in details 
the characteristics of the source during high and low count 
rate, we separated the e class light curve in two distinct 
types where type I is characterized by higher count rate and 
type II has lower count rate. We also study the hardness- 
intensity diagram by studying count rate in three different 
energy bands. In hardness-intensity diagram we have plotted 
HR2 vs Total flux {X + Y + Z) in three bands. To investi- 
gate any presence of low-frequency quasi-periodic frequency 
(LQPO) between 0.1-10 Hz, we produce power-density spec- 
trum by using powspec. The power spectrum is normalised 
such that their integral gives the squared rms fractional vari- 
ability and the expected white noise level is also subtracted. 
We use 10 milli sec. time-resolved data with geometrical bin- 
ning of -1.06 for the study of power-density spectrum. To 
study variation of rms amplitude, we extracted light curves 
between 2.0-8.0 keV and 16.0-60.0 keV using binned data 
having 0.1 sec time resolution. From these, we have calcu- 
lated RMS amplitude for these two different energy bands 
for each 1 sec binning. 

We have extracted spectrum from PCA as well as 
HEXTE data. For PCA data reduction, we have used all 
the cleaned data file having 16s resolution and correspond- 
ing response matrix file has been created adding all xenon 
layers of PCU 0, 1 and 2 using pcarsp. Then from the same 
data file we have extracted background spectrum using run- 
pcabackest. For HEXTE data reduction, we have extracted 
source and background file using hxtback. Dead time cor- 



3 THE NEW CLASS e 

In the upper panel of Figure 1, RXTE All Sky Monitor 
(ASM) light curve, at the time of the occurrence of this 
new class is plotted and in the lower panel of same Figure 
Green Bank Interferometer (GBI) radio flux density of the 
source is plotted for the same time. In both of these flgures 
the time of new class is marked in blue color. It is found 
that the new class took place when the source was in radio 
quiet state, weeks before a radio flare. The source was in a 
X-ray flare for all the time when the class was observed and 
always it was observed in the rising part of the flare. 

3.1 Timing features of the new class 

The new class is characterized by nearly periodic fluctua- 
tion between higher and lower count rate where the higher 
count rate is about 3-4 times more than the lower count rate. 
In the upper panel of the Figured a typical light curve of 
GRS 19154-105 at the time of e class is shown, which was 
observed on 23rd April 1999. The source shows higher count 
rate for 350 to 500 s followed by lower count rate of 80 to 100 
s. This indicates some kind of thermal-viscous instability in 
the disk which is adjusting itself to accommodate higher ac- 
cretion rate at the time of X-ray flare. In the lower panel of 
the Figure[21 variation of HR2 with time is shown. It is found 
that at the time of type II, HR2 increases to ~ 0.5 from ~ 0.1 
of type I and thus while lightcurve is varying between type 
I and type II, the source is undergoing an oscillatory state 
transition between soft and hard states. In the upper-panel 
of Figure [3] we plotted color-intensity diagram, where the 
variation of total flux with HR2 is shown. In the lower panel 
of Figure [S] we have shown color-color diagram, where vari- 
ation of HRl with HR2 is shown. In color-color diagram and 
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Figure 2. Upper panel: A typical lightcurve of GRS 1915+105 
in newly discovered e class in the energy range 2-60 keV. The 
observation took place on 23rd April 1999. Lower panel: Variation 
of HR2 for the same time as in upper panel. 



hardness-intensity diagram, points corresponding to type I 
and type II are shown in green and red color respectively. 
In these diagram type I and type II form two clearly sepa- 
rable regions. Since the intensity in type I is high and HR2 
value is low, it is in high-soft state. The type II region is in 
the hard-intermediate state. In Figure 4 , the power-density 
spectrum of the source is shown. No, definite low-frequency 
quasi-periodic oscillation (LQPO) has been observed. After 
fitting with lorenzian between the frequency range of 0.01 to 
10 Hz, we calculated total power 0.016 for type / and 0.116 
for type II. Total power of type I state is ~ 7 times lower 
than the total power of type II, this result is also expected if 
we assume type I to be soft state and type II to be hard state. 
From our analysis it is clear that the type II state is mostly 
similar to state 'C while the type / state is mostly similar to 
sta te 'B' where d efinition of the staes are same as mentioned 
by iBelloni et al. I (|2000. ). So, it is clear that transition from 
state 'C to state 'B' and and state 'B' to state 'C both are 
possible. This class is special as in the previous study with 
all earlier 12 classes iBelloni et al.l (|2000l ) demanded that the 
transition from state 'C to state 'B' is not possible and e is 
the only class, so far discovered, where this kind of transition 
is possible. 

To check the variability of different ban ds in e class, we 
plot RMS amplitude variation in Figure 5 . iKalemcil (|2002| ) 
has showed that rms amplitude of variability increases with 
energy for a variability system with significant disk emission 
but the relation is insignificant when disk emission is suffi- 
ciently low. Our RMS amplitude diagram clearly show that 
average RMS variability for e class is significantly high for 
high-energy band (14.0-60.0 keV) compared to low energy- 
band (2.0-8.0 keV) in case of type-I. From spectral analysis 
we know that type-I spectra is mostly dominated by disk 
blackbody component. Hence for e class, the proposition 
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Figure 3. Upper panel: hardness-intensity diagram, showing 
variation of total flux with HR2 for the data observed on 23rd 
April 1999. Lower panel: color-color diagram, showing variation 
of HRl with HR2. In both of these figures, points corresponding 
to type I is plotted in green and points corresponding to type II 
is plotted in red. Points corresponding to type I and type II form 
two clearly distinguishable separate region on both of the plots. 
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Figure 4. Power density spectrum for type / and type //observed 
on 23rd April 1999. Total power of type II state is higher than 
that of type I. 

that presence of black body component in spectra reduces 
the rms variability in lower energy band is clearly seen. 

3.2 Spectral features of the new class 

The spectral fitting parameters for type I and type II is 
shown in Table 1 and 2. The 90% confidence level for each 
fitted parameter and x^/f s^re also shown. The flux of the 
individual components are also included. In Figure 6 , we 
have shown different components for type I and type II of 
the spectrum. 

From the best fitted values, it is found that type / spec- 
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Table 1. Observation ID and Spectrum fitting parameter for Type I 



Obs. ID Date Ti r2 Break x^/f F2-26 Fpow Fblackbody Firo: 

Energy (keV) (koV) 



40403-01-07-00 


04/23/1999 


3.42±0.06 


4.81±0.06 


18.3±0.3 


1.95±0.01 


1.03 


2.16 


0.81 


1.35 


0.004 


40703-01-13-01 


05/01/1999 


3.58±0.07 


4.87±0.06 


17.5±0.3 


1.99±0.02 


0.96 


3.22 


1.64 


1.57 


0.006 


40703-01-14-00 


05/08/1999 


3.21±0.04 


4.77±0.06 


16.7±0.3 


1.90±0.00 


0.92 


2.77 


1.25 


1.52 


0.004 


40703-01-27-01 


08/23/1999 


3.11±0.06 


4.91±0.10 


16.0±0.3 


1.94±0.01 


0.96 


1.82 


0.74 


1.09 


0.001 


40703-01-27-00 


08/23/1999 


3.54±0.05 


4.92±0.11 


18.1±0.4 


1.97±0.01 


1.03 


2.41 


1.07 


1.33 


0.006 


40703-01-28-02 


08/28/1999 


3.33±0.06 


4.88±0.11 


19.0±0.4 


2.01±0.02 


0.85 


1.71 


0.59 


1.12 


0.004 


40703-01-29-00 


09/05/1999 


3.28±0.04 


4.77±0.06 


17.8±0.3 


1.92±0.01 


1.06 


2.66 


1.25 


1.40 


0.005 


Column density 


^ 7 X 10^2 cm" 


^ was used. Flux is given 


in the unit of photon (1 photon 


= 9.012 X 10^8 


ergs /cm" 


"^/sec in 


2-26 keV 



band). 

Table 2. Observation ID and Spectrum fitting parameter for Type II 



Obs. ID Date Tl Tt X^/f F2-26 Fpow Fuackbody Firo, 

(keV) 



40403-01-07-00 


04/23/1999 


3.58±0.03 


2.06±0.02 


0.98 


1.15 


0.72 


0.43 


0.002 


40703-01-13-01 


05/01/1999 


3.65±0.04 


1.97±0.02 


0.96 


1.94 


1.26 


0.67 


0.008 


40703-01-14-00 


05/08/1999 


3.63±0.04 


1.99±0.02 


1.00 


1.09 


0.64 


0.44 


0.005 


40703-01-27-01 


08/23/1999 


3.64±0.04 


2.02±0.02 


1.06 


1.06 


0.60 


0.47 


0.004 


40703-01-27-00 


08/23/1999 


3.71±0.05 


2.00±0.02 


1.11 


1.10 


0.62 


0.47 


0.008 


40703-01-28-02 


08/28/1999 


3.54±0.04 


2.07±0.18 


0.98 


1.18 


0.61 


0.57 


0.003 


40703-01-29-00 


09/05/1999 


3.52±0.03 


1.98±0.02 


0.94 


1.32 


0.69 


0.62 


0.004 



Column density ~ 7 X 10^^ cm ^ was used. Flux is given in the unit of photon (1 photon= 9.012 X 10 ergs/cm ^/sec in 2-26 keV 
band). 
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Figure 5. Top Panel: Variation of RMS amplitude in e class for 
16-60 keV band as observed on 23rd April 1999. middle panel: 
Variation of RMS amplitude for 2-8 keV band. In the bottom 
panel we show ratio of the RMS amplitude for the above two 
mentioned band. 



only 26% to 43%. Hence, it is cl ear that it is quite similar t o 
HS (High State) as described bv lHoman and I3ellonil l|2005h . 
The break in the power-law continuum suggests that there 
is two-temperature accretion flow in this mode. In case of 
type II spectra, it is found that the hard power law com- 
ponent is contributing 58% to 66% of the total flux, while 
the disk black body is contributing from 40 to 43% of the 
total flux. Hen ce type I is similar t o the LS (Low State) 
as described by iHomaJi and BeUonil (|2005l ). It is a very im- 
portant fact to notice the temperature of diskbb is very high 
and normalization value is low, suggesting very hot disk and 
very low inner radius. Hence, accretion rate is very high. We 
have also tried to fit the spectrum with Comptt model, but 
the contribution of compton cloud flux to the total flux is 
less than 1% and the optical depth is very low, suggesting 
the absence of optically thin cloud. Looking at the spectrum 
point of view, we can conclude that this e class is basically 
intermediate transient type of class where the power-law and 
thermal component have varing relative contributions. 



tra can be fitted well with broken power law with a break 
near 18 keV, along with disk black-body model and small 
Gaussian at ~ 6.4 keV while all the type 77 spectra are found 
to be fitted with simple power-law, along with disk black 
body and small Gaussian. The contribution for individual 
modeling components are also shown in the same figure. For 
type I spectra, the contribution of thermal component fiux 
to the total fiux is varying between 56% to 70% between 
different observations, while the power-law fiux contributed 



4 DISCUSSION 

The class is observed when GRS 1915+105 was undergoing 
rising phase of a X-ray fiare though the source was quiet 
in radio. Studying the new class sequentially, it can be con- 
cluded that as the X-ray flare intensity is getting higher, the 
e class is spending longer time in type I and lesser time in 
type 77 and gradually the type 77 vanishes. The phenomenon 
inevitably suggested that the e class is an oscillatory transi- 
tive stage between low and high state in the color-intensity 
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Figure 6. PCA spectrum of GRS 1915+105 in e class observed on 23rd April 1999. In left panel the spectrum and different model 
component in type I region and in right panel spectrum and different model component in type II region is shown. 



evolution diagram and finally get stabilized at hard state. 
This show that at time of e class, the source was undergoing 
some thermal viscous instability when it was trying to ad- 
just between hard and soft states when the strong emission 
from radio jet was not ON and the source is in fully radiative 
mode. 



This paper has been typeset from a T[5X/ WV^^ file prepared 
by the author. 
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